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Abstract 



^ ■ We consider the total production cross section of heavy coloured particle pairs in 

hadronic collisions at the production threshold. We construct a basis in colour 
space that diagonalizes to all orders in perturbation theory the soft function, 
which appears in a new factorization formula for the combined resummation of 
soft gluon and Coulomb gluon effects. This extends recent results on the structure 
of soft anomalous dimensions and allows us to determine an analytic expression for 
the two-loop soft anomalous dimension at threshold for all production processes 
of interest. 



1 Introduction 



Perturbative calculations of partonic cross sections at hadron colliders often fail near the 
boundaries of partonic phase space due to logarithmically enhanced terms from soft gluon 
radiation. If it can be argued that the hadronic cross section is dominated numerically by 
these threshold logarithms, they should be summed to all orders in perturbation theory. 
This can be done in Mellin moment space [1-5] or directly in momentum space [6-8]. In 
either case, the theoretical basis for resummation is a factorization of the partonic hard- 
scattering cross section a in the partonic threshold region into hard and soft contributions 
of the schematic form 

a = H(g)S (1.1) 

with a hard function H and a soft function S. The soft function satisfies an evolution 
equation, whose driving term is the anomalous dimension matrix of S. Resummation 
amounts to solving this equation. 

Particularly important for experiments at LHC and Tevatron are pair production pro- 
cesses of heavy coloured particles H, H' in a collision of hadrons N± and N 2 , 

N 1 (P 1 )N 2 (P 2 ) -> H( Pl )H'(p 2 ) + X. (1.2) 

In this case the partonic cross section contains terms of the form [a™ ln m (3 ] ( "threshold 
logarithms") and (a s / ' (3) n ("Coulomb singularity"), where (3 = (1 — 4M 2 /s) 1 / 2 is the heavy 
particle velocity, which are enhanced near the partonic threshold s ~ 4M 2 , with M the 
average heavy-particle mass. The threshold logarithms have been discussed in the past for 
various production processes of the form (jl.2p leading to improved predictions for the top- 
quark production cross section at hadron colliders [9-17], production of supersymmetric 
coloured particles [18-20], and colour-octet scalars [21]. In two-to-two scattering the soft 
function is a matrix in colour space. In general this matrix depends on the kinematical 
invariants of the scattering process [3]; however, for resummation of threshold logarithms in 
the total cross section, the relevant quantity is the soft anomalous dimension at threshold. 
In this case, the question arises whether the resummation of logarithms is altered by the 
presence of Coulomb corrections, which must also be summed. This has not been addressed 
in the past, where the Coulomb correction is technically considered as part of the hard 
function H. 

In a separate paper [22] we employ effective field theory and field redefinitions to derive 
an extended factorization formula of the form a = H <g> J <g> S for the hard-scattering 
total cross sections related to (II. 2p . which implies a proof of factorization of soft gluons 
in the presence of Coulomb exchange. More precisely, the cross sections of the partonic 
subprocesses p(ki)p'(k 2 ) — > H(pi)H'(p 2 ) + X where pp' G {qq,qq, gg, gq, gq} are expressed 
as 

o pp ,{s,n) = Y J H ll >{M^) J du JR a (E-^)W^(uj^). (1.3) 

Here s is the partonic centre-of-mass energy, and E = \/I — 2M. The formula applies 
to heavy particle pairs produced in an S-wave and is valid at the leading order in the 
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non-relativistic expansion. The new function J# Q sums Coulomb gluon exchange related 
to the attractive or repulsive Coulomb force in the irreducible colour representations R a 
that appear in the product representation R (g) R' of the final state particles and includes 
the leading Coulomb singularities (a s //3) n . We note that Coulomb summation has been 
included in various forms in [16,17,19] where, however, the factorization of Coulomb from 
soft gluons is put in as an assumption. A formula equivalent to eq. (jl.3p for the factorization 
of electromagnetic effects in P^-pair production at e~e + colliders has been derived in [23]. 
The purpose of the present paper is to discuss the colour decomposition of the generalized 
soft function W Ra , which is the crucial ingredient for resummation, and to provide the two- 
loop anomalous dimensions and one-loop soft functions, which are necessary for next-to- 
next-to-leading logarithmic (NNLL) resummations. We shall also briefly outline the general 
structure of resummation and explain further NNLL effects not included in eq. (11.31) . 

The soft functions are vacuum expectation values of Wilson line operators that retain 
only information about colour and direction of massive and light-like particles in the hard 
process. The light-like (n 2 = 0) Wilson line for an incoming particle in the representation 
r of SU(3) with generator T (r)a is 



S%)( x ) = Pexp 



Ws 



o 



dsn- A a (x + ns)T {r)a 



:i-4) 



For an outgoing heavy particle in representation R we define the time-like (v 2 = 1) Wilson 
line 



Si R)t (x) = Pexp 



'(]> 



poo 

/ dsv A a {x + vs)T {R)a 
Jo 
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The inverse (adjoint) Wilson line operators follow from replacing path-ordering by anti 
path-ordering, and ig s — ► — ig s . The soft function W Ra (uj,n) in (11.31) is a descendant of a 
more general soft function defined by 

W$ } (z,ii) = (Q\T[S v 

Mk2Svtek 1 Slj b2 S^ ib ^(z)T[S ntaii S n ,a2jSi,k 3 a 3 S^^ (1-6) 

where T and T denote time-ordering and anti time-ordering, respectively, n and n denote 
two light-like vectors satisfying n ■ n = 2. The superscript on the Wilson lines denoting 
the colour representation has been omitted. In the factorization formula (11.31) we need the 
Fourier transform of the soft function defined according to 

WiSMA = J d ^e^ 2 W^ b \(z Ari. (1.7) 

In (11.31) we also decomposed the colour multi-indices in a set of basis structures cf\ defined 
below, and performed a projection on the irreducible representations R Q of the final state 
particle pair. Thus W Ra (u,{i) is given by 

WS-^/i) = ^Hg( W)/ i) C g*. (1.8) 
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The remainder of the paper is concerned with the properties of the soft functions (11.61) and 
(11. 8p . We remark that the soft function that appears in (11. ip follows from contracting (11. 6p 
with the trivial colour factor 5k 1 k 3 Sk 2 k 4 instead of the projectors on the representations 
of the heavy-particle pair as in (11.81) . Therefore the factorization of the Coulomb gluons 
results in a more complicated colour structure of the soft function than the original one, 
Wu>, which is simply W?, a due to the completeness of the projectors. For this reason, 
the original factorization formula (ll.ip should only be used for partonic thresholds where 
the relative velocities of the final state particles are relativistic, and Coulomb exchange 
is not enhanced. In this situation, the relevant soft function is (11.61) with two unequal 
four- velocity vectors Vi, v 2 for the final state particles, contracted with 5k 1 k 3 8k 2 k 4 ,- The 
main simplifications discussed in the following hold only for the soft function at threshold 
(t>i = f 2 ), and are thus applicable to the total partonic cross sections. 

The paper is organized as follows. In section [2] we discuss the colour structure of the 
soft function. The simple kinematical structure of the soft function W Ra at threshold 
allows us to construct a colour basis that diagonalizes the soft function to all orders of 
perturbation theory. We shall see that the four-particle soft function can be reduced to 
three-particle soft functions corresponding to a single heavy particle in the final state. In 
section [3j after discussing the general structure of soft-gluon and Coulomb resummation 
near threshold, we calculate the one-loop soft function for general representations of initial 
and final state particles and obtain the two-loop soft anomalous dimensions employing 
results from [25]. This supplies all ingredients for NNLL resummations for top-quark or 
sparticle pair production associated with the (leading) soft function (11.81) . The two-loop 
soft anomalous dimension at threshold exhibits Casimir scaling as has been found explicitly 
at one-loop in previous examples [3,10,18]. Technical results related to colour, Wilson lines 
and anomalous dimensions are summarized in a number of appendices, including the colour 
bases for triplet, anti-triplet and adjoint coloured particles in the initial and final state, 
which covers the cases of interest. 

2 All-order colour structure of the soft function 

The physical picture of production of a heavy particle pair at threshold suggests that soft 
gluon radiation cannot resolve the two particles and couples to the total colour charge of 
the pair, determined by the representation R a in the decomposition of the product repre- 
sentation R ® R'. Therefore the structure of the (leading) soft function for the production 
of a non-relativistic particle pair in a given representation should be that of a single par- 
ticle in the same representation. This is in agreement with the result that the one-loop 
soft anomalous dimension at threshold is proportional to the quadratic Casimir operator 
of the representation of the heavy-particle pair [3, 10, 18]. In this section we show how to 
obtain this structure, to all orders in perturbation theory, from the expression (ll.8p for the 
soft function. In order to accomplish this, we construct the projection operators P Ra and 
the basis tensors c- 1 ' in the definition of the soft function (11. 8p from Clebsch-Gordan coef- 
ficients for the decompositions of the product representations of the initial and final state 
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systems into irreducible ones. In this basis it is then easy to show that the soft function is 
diagonal to all orders of perturbation theory. Some technical details of this construction 
are relegated to appendix [A] Explicit colour bases and projectors for all production pro- 
cesses of heavy particles in the fundamental, antifundamental and adjoint representations 
are provided in appendix [B] This covers all production processes of squarks and gluinos 
and completes previous results for the colour bases of squark-antisquark and gluino pair 
production [18]. 

2.1 Notation 

With respect to colour, we use a notation that does not distinguish particles and antipar- 
ticles. If, for example, H' is the anti-particle of H, then this convention implies that H' 
transforms in the complex conjugate SU(3) representation of H . Similarly, an initial-state 
antiquark transforms in the anti-fundamental representation. The Wilson lines S inherit 
the corresponding representations. It is useful to perform a decomposition of the product 
of the representations of the final state and initial state particles into irreducible represen- 
tations: 

r©r' = ^r Q , R®R' = ^R a . (2.1) 

a R a 

Examples relevant for the standard model and most extensions of the standard model are 
particles in the fundamental and the adjoint representation where we have the decompo- 
sitions 3 © 3 = 1 + 8 (e.g. top-antitop and squark-antisquark production), 3 © 3 = 3 + 6 
(e.g. squark-squark production), 3©8 = 3 + 6 + 15 (e.g. squark-gluino production) and 
8®8 = 1©8 S ©8 A ©10©I0©27 (gluino pair production). Generators of the SU(3) trans- 
formation in a representation R are denoted by T^ a . In practice we need the generators 
in the fundamental and anti-fundamental representations, T^a 2 = T£ , = — T® , 

and in the adjoint, T£fi£ = if aiaa2 = F£ a . The quadratic Casimir operator for a repre- 
sentation R is denoted by (T^ a T^ a ) aia2 = C R S aia2 . 

2.2 Construction of projectors on irreducible representations 

We now review a general method to construct the projection operators on the irreducible 
representations of the heavy-particle pair appearing in the soft function (11.81) from the 
Clebsch-Gordan coefficients that combine two objects transforming in the representations 
R and R' of the group into a single object in an irreducible representation R a . The 
Clebsch-Gordan coefficients are defined as a unitary basis transformation, 

e aia 2 = ^ ] ^aa 1 a^ e a"; (2-2) 
Ra 

from the basis vectors of the tensor product space R®R' to basis vectors of the irreducible 
representations R a (see e.g. [26]). Here the e ai(l2 are elements of a complex, orthonormal 
basis of the tensor product R © R! while the e^ a are elements of a basis of the irreducible 
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representation R a . In this equation and in the following, repeated indices are summed 
over. Unitarity of the basis transformation implies 

Cg£* a2 bo = S aibl 5 a2b2 , (2.3) 



Ra 



Caa*a2 C ^a ia2 ~ $R a R $al3- (2.4) 



A vector V in the tensor-product space can be written in the two bases as V = V aia2 e aia2 = 
J2n a Va e a, where the components are related by 

V ai a 2 = J2 C ^: a2 V a . (2.6) 

Ra 

( Fi) 

Consistency with the group transformations V a — > U^JVp and the corresponding transfor- 
mations of V aia2 implies that the Clebsch-Gordan coefficients satisfy 

r*R a tt( r ) n( R,S) — tt( r <*) r^Ra (n j\ 

, ^aa 1 a 2 u a 1 b 1 u a 2 b 2 ~ u af3 °/36i& 2 - V '/ 

We can now construct the projectors P Ra on the irreducible representations from the 
coefficients C Ra according to 

a\a 2 ao,a^ ^ aaia 2 ^~ / aa^a^' V / 

The orthonormalization condition (12. 4p implies that these are projectors satisfying 

pRa p R P _ X pRa (2 Q) 

a\a 2 b\b 2 b\b 2 c\c 2 u R a Rp a\a 2 c\c 2 i 

and the relation (12.31) is equivalent to the completeness relation 

yi^hh = ^161^262- ( 2 - 10 ) 

Ra 

As an example, consider final state particles in the 3 and the 3 representation of SU(3), 
e.g. top-antitop and squark-antisquark pairs. The two particles combine to a singlet or 
octet, and the Clebsch-Gordan coefficients read 

C£ 2 = -4. Ka 2 , Cffi « = V2T« ai . (2.11) 



We obtain the familiar projectors 

pi 1 ) _ J_ X X p(8) _ .)■.». 7 ' (•'>'{•'>) 
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2.3 Colour basis for the hard production process 

The colour structures $K in the definition (ll.8p describe the colour structure of the hard- 
scattering process. They provide a decomposition of the scattering amplitudes of the 
partonic process pp' — ► HH' into independent basis tensors: 

-4pp'{a} = 5Z C {a}*^PP" (2.13) 

i 

Here we use a multi-index notation {a} Repeated multi-indices are summed 

over all four components. We take the cW as an orthonormal basis satisfying 

( 2 - 14 ) 

The colour basis structures can be chosen as invariant tensors in the representation r <S> 
r' ® -R ® -R' - 

_(») _ rrtflOtr/WtJO r/W rrM fnic^ 
C {a} — ^363^464 C {fe} U b iai U b 2 a 2 - \*- LO ) 

A colour basis that is convenient for the discussion of the colour structure of the soft 
function can be constructed from the Clebsch-Gordan coefficients, similarly to the projec- 
tors. Consider the subset of representations that appears in both the sets {r a } and {R a } 
in (12. ip . treating multiple occurrences of equivalent representations in the decomposition as 
distinct. We then form pairs Pi = (r a , Rp) of equivalent representations r a and Rp, where 
the index i enumerates the allowed combinations. For example, in case of8®8^8®8 
we have the eight pairs 

Pi e {(1,1), (85,85), (8 A , 8 S ), (8 A ,8 A ), (8 S ,8 A ), (10,10), (T0,10), (27,27)}. (2.16) 

The 10 and 10 are inequivalent representations, so the pair (10, 10) is not allowed. For the 
allowed pairs Pi, the colour structures 

„(*) _ \ nr a ry R /3* Icy i ? \ 

y/<hm(r a ) 

form an orthonormal basis satisfying (12.141) . By construction, the operators are invariant 
tensors under global SU(3) transformations satisfying (I2.15p . In appendix |X] we use colour 
conservation of the scattering amplitude to show that the basis tensors always can be 
chosen as in f)2.17p . 

To illustrate this construction, we consider heavy particles in the 3 and 3 produced 
from a quark- ant iquark initial state. There are only two possibilities to combine the initial 
and final state representations, Pi = {(1, 1), (8,8)}. Using (12. lip we obtain two operators 
for the basis of the hard-scattering process: 

M _ _/y U 0l £l 2 L '0304) c {a} ~ fjy- J -a 2 a 1 J -a :i a 4 l z - io J 



with Da = N% — 1. This is the same basis found to diagonalize the one- loop soft anomalous 
dimension at threshold [3]. A complete list of all projectors and basis tensors for the 
production of particles in the fundamental and adjoint is given in appendix IBl 
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2.4 Diagonalization of the soft function 

We now show that the soft function is diagonal to all orders of perturbation theory in 
the basis constructed in section 12.31 To achieve this, we express the components of the 
soft function (11. 8ft in terms of Wilson lines in the representations R a , as suggested by the 
physical picture of soft radiation off the total colour charge of the final state system. We 
first note that the Wilson lines satisfy a relation analogous to (12.71) . 

ryRa q(R) q( R ') _ q( R ct)s-iR a /q iq\ 

^aa 1 a 2 °v,a 1 b 1 °v,a 2 b 2 ~ °v,af3 ( - y f3b 1 b 2 - l z - ±y J 

A proof of this relation is given in appendix El Using this identity and the completeness 
relation ( 12.31) of the Clebsch-Gordan coefficients we can write the soft function (11.61) in 
terms of Wilson lines for single particles in the irreducible representations R a : 

wWfV ,,\ _ n Rf} * n Rf) r Ra r Ra * 

vv {ab}\ z i AV — /^i °/96364 /cfeifc2^aa 3 a4 Ly A,A;3fe4 

Ra,Rp 

x (0\T[SX S UA .ttxlWT^^^KOJIO). (2.20) 

This combination of the two Wilson lines in the representations R and R' to a single one 
in the representations R a is only possible if the heavy particle pair is produced close to 
threshold and not for generic kinematics where Wilson lines in different directions Vi,v 2 
appear in the soft function. 

In (12.201) the Wilson lines related to the final state system are still in two different 
representations. This structure simplifies when we compute the components of the soft 
function (11.81) where the soft function is contracted with a projector on an irreducible 
representation. Using the definition of the projectors (12. 8p we have the identity 

C KMk 2 P {k} C \teki = A &R a i?, 7 $R a Rp j (2.21) 

i.e. the projectors enforce that the representations Rp and R 1 are identical. For distinct 
equivalent representations such as the 85 and 8^ representations in the decomposition 
of 8 ® 8 we obtain a vanishing result. We can then express the components of the soft 
function (jl.8p as 

W«r(u,,ri = c^Wfc^^p*, (2.22) 

where we have introduced the soft function for the production of a single particle in the 
representation R a 

ftfemM = (0im^4 j62 <*J(^)T[^, ail ^,a 2 ^S](0)|0). (2.23) 

Here we have extended our multi-index convention to the indices of the irreducible repre- 
sentations of the final state system by defining {aa} = a\a 2 a. Note that the index k of the 
two final state Wilson lines in (12.231) is contracted, so this is analogous to the soft function 
appearing in the conventional treatment where Coulomb gluons are not factorized. This 
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Figure 1: Structure of the simplified expressions of the soft function (12.221) and (I2.23p . 



function has been considered in [21] for the case of a single colour-octet scalar. The colour 
basis tensors of the production operators in the new notation are given by 

R a (i) _ (i) r R a _ 1 nr a r (n ja} 

°{aa} — C {a}^aa 3 a 4 — rp 7 rVa 2 " fi « R ?' K^-^l 

11 11 A/dim(r a ) 

As indicated, they are nonvanishing only if the final state representation R a is identical to 
the final state representation Rp in the pair Pj = (r a , Rp) that defines the tensor This 
implies that r a and R a must be equivalent, as is intuitively clear. 

In ( 12.23P we have reduced the problem of soft gluon radiation in pair production of 
heavy particles to that of the production of a single particle in the representation R Q . This 
is the first main result of this paper and will simplify the computation of the one-loop soft 
function and the two-loop soft anomalous dimension. The structure of the results (I2.22p 
and (I2.23P is sketched in figure [TJ 

Eq. (I2.24p for the basis tensors automatically leads to the vanishing of a large number 
of components of the soft function. The component W Ra is nonzero only if the repre- 
sentation R a of the final state system is identical to both final state representations Rp 
and Rp and hence, by construction of the basis tensors, equivalent to both initial state 
representations r a and r a > in the pairs Pj = (r a , Rp) and p/ = (r a /,Rp/). Therefore the 
soft function is block-diagonal in the basis (I2.17P with off-diagonal elements arising only if 
several representations in the decomposition (12.1 p of the initial state system are equivalent. 
(A non-trivial structure does not appear if the several equivalent representations appear 
only in the final state since the projectors P Ra always project on a unique representa- 
tion.) Since in the decompositions of 3 ® 3, 3 <S> 3 and 3 ® 8 no representation occurs more 
than once, the matrices W Ra are diagonal if at least one quark or anti-quark is present in 
the initial state. The only example where a non-trivial matrix structure can arise is the 
gluon- gluon channel. 

For the example of an 8 <E> 8 final state produced from gluon-gluon fusion, according 
to (I2.16p . only the soft functions for the two octet final states 85 and 8,4 are (potentially) 
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non-trivial two-by-two matrices, 



W 8s 



/o 



\ 



/o 



w 



8,1 



■J 



\ 



wit wl- 
wit w*£ 



\ 



(2.25) 
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with entries in the order of the basis elements used in (12. 16p . The soft functions for 
the 1, 10, 10 and 27 representations, however, consist of a single non-vanishing element. 
Analogously, for the production of a 3 <S> 3 final state from gluon-gluon fusion the singlet 
soft function consists of a single non-vanishing entry while the colour-octet soft function 
contains a potentially non-trivial two-by-two submatrix mixing 8a and 85 initial states. 

We now use Bose symmetry of the soft function to show that the off-diagonal elements 
in the soft function for the production of a colour-octet state vanish, so that the matri- 
ces (12.251) and the analogous ones for a 3 <S> 3 final state are in fact diagonal to all orders in 
perturbation theory. First note that, independent of the nature of the final state system, 
these off-diagonal matrix elements involve a combination of Wilson lines contracted with 
the Clebsch-Gordan coefficients Caaia 2 and Caaia 2 defined in ( IB. 51) : 

^/3fri6 2 [^n,ih ^n,jb 2 }( Z ) [Sn,a 2 jSn, ai i\ (0)C,^1 2 OC [S^^D^^ SljfJ (z) [Sn,a 2 jF" 2ai Sn, ai i\ (0) . 

(2.26) 

Because the two incoming Wilson lines are indistinguishable, the soft function must be 
invariant under the exchange of their colour labels and momenta, i.e. the exchange (n, i) <-> 
(n,j). This statement translates into the following equation: 

[^n,ibi ^6i6 2 Sftjhi] ( 2 ) [^n,a 2 jF a2ai S n , ai i] (0) 

= [^n,ib 1 F > b 2 b 1 ^n,jb 2 \( Z ) [^n,a 2 jF aia2 ^n,aii\ (0) 
= — [^ri,i&i A>ib 2 ^n,i*J ( 2! ) \Sn,a 2 jFa 2ai S nAl i\ (0) 

= 0. (2.27) 

Here we have used the symmetry properties of the D and F tensors. Therefore the off- 
diagonal terms of the colour-octet soft function such as in (12.251) vanish. Furthermore from 
the explicit results in section [3] we deduce that the one- loop soft function and two- loop soft 
anomalous dimensions are determined by the quadratic Casimir operators, so the diagonal 
elements of W 8a and W 8s are all identical. This would extend to higher orders if Casimir 
scaling held to all orders. However, presently we do not have a proof for this. 

To summarize we have shown that for all initial states relevant to hadron colliders the 
soft function is diagonal to all orders in a colour basis for the hard-scattering process given 
by (I2.17p . This holds independent of the nature of the final state particles, i.e. equally for 
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top-quarks, squarks or gluinos. As mentioned at the end of section [TJ the approach used in 
most phenomenological applications corresponds to using the soft function Wa> obtained 
by summing up all final state representations so our results apply to this case as well. Let 
us briefly recall the main ingredients used in order to arrive at this result: 

• The Coulomb interaction is diagonalized by the decomposition of the final state sys- 
tem into irreducible representations (12. ip . leading to the definition of the components 
of the soft function (11.81) where the Wilson lines associated to the final-state particles 
are projected onto the irreducible representations. 

• The hard-scattering amplitudes are colour conserving, eq. ( 12.1 5ft . which allowed us 
to choose the basis of colour tensors according to (12. 1 TP . 

• For a heavy particle pair produced directly at threshold both particles have the same 
velocity, allowing to combine the two final-state Wilson lines into a single one. 

• Due to Bose symmetry of the soft function there is no interference of the production 
from a symmetric and antisymmetric colour octet. 

2.5 Examples 

In this subsection we give some examples of the formalism in order to show how it is related 
to the colour bases used in previous computations. As an illustration of the result ( I2.22j) 
for the components of the soft function, we give explicit expressions for the example of a 
3 <S> 3 final state and a quark-antiquark or gluon-gluon initial state. In subsection 12.5.21 
we compare the basis for an 8 ® 8 final state produced in gluon fusion to that used previ- 
ously [18]. 

2.5.1 Soft function for a 3 ® 3 final state 

For quark-antiquark initiated production of a heavy particle pair in the 3® 3 representation, 
the required colour basis is given by (12. 18ft so the soft functions for the singlet and octet 
final states, W^, and Wf v , are two-by-two matrices. The only non- vanishing component of 
the soft function for the singlet case is identical to that in Drell-Yan production [8, 27] 

Wl^n) = ^(OlTrfT^^T^^^liO) = W DY (z^). (2.28) 

Here we have deviated from our usual notation of using the anti-fundamental representation 
for anti-particles and expressed the result using only Wilson lines in the fundamental 
representation in order to simplify the matrix structure. Similarly, the only non-vanishing 
component of the colour-octet soft function is given by 

W» 2 (z,pi) = -j^-^ (0|Tr[T[^tr^](^)T[4T^ n ^J(0)]|0), (2.29) 
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where the S v are in the adjoint representation and the trace is over the fundamental 
representation. 

For the production of a 3 <S> 3 final state from gluon fusion there are three possible 
combinations of initial and final state representations: 

Pi e {(1,1), (8 5 ,8), (8 A , 8)}. (2.30) 

The Clebsch-Gordan coefficients and colour-basis elements for this case are collected in 
appendix ( IB. II) . Since the set of Pj's has three elements, the singlet and octet soft functions 
are three-by-three matrices. The only non-vanishing element of the soft function for the 
singlet channel is, up to normalization, again given by (12.281) . where now the Wilson lines 
S n and S n are in the adjoint representation. An octet final state can be produced either 
from a symmetric or antisymmetric octet initial state corresponding to the basis elements 
c*- 2 -* and in ( IB. 61) . The non- vanishing diagonal elements of the soft-function matrix for 
the octet channel in this basis are given by 

WUz^) = — 2 _ ^ _ (0|Tr[T[5' UjOC S'tD a 5 n ](^T[4DX< c6 ](0)]i0), (2.31) 

WU^V) = AT( 1 T\ (0|Tr[T[^ ifec StF 6 ^](^)T[4P a ^St iC J(0)]|0). (2.32) 
W - 1) 

Here we have used the fact that the Wilson lines in the adjoint, S n , S n , are real to write 
the matrix product in a convenient form. The expressions for the octet soft functions 
agree precisely with [21] (up to their notation for the result in momentum space). The 
off-diagonal elements and involve the structure (I2.26P and therefore vanish by 
symmetry arguments. 



2.5.2 Colour octet states in8<g>8-H>8<g>8 

We would like to comment briefly on previous results for the basis for the gluon-induced 
production of heavy particles in the adjoint representation (e.g. gluinos). The complete 
basis for this case is given in appendix [B] Here we will only need the four operators 
corresponding to the different combinations of 83 and 8^ in the initial and final state (IB. 201) . 
This basis differs slightly from the one constructed in [4] for dijet production where the 
linear combinations (c^ ± c' 5 ') have been used. In [18] that basis has been shown 
to diagonalize the one-loop soft anomalous- dimension matrix corresponding to the soft 
function Wn> discussed below (11.81) at threshold, i.e. for v\ = v 2 - From the general 
arguments given above and from an explicit calculation we find that the one-loop soft 
functions for colour octet final states W~f and W»f are diagonal in the basis (IB. 201) but 
not in the one used in [18]. However, since the Coulomb functions J 8a and J 8s are identical, 
only the sum W Ss + W 8a enters the cross section (II. 3p . Since the off-diagonal terms cancel 
in the sum, our result is consistent with the one in [18]. Similar remarks apply to the 
colour tensors related to the 10 and 10 where the basis used in [4, 18] is appropriate for the 
sum W 10 + W 10 that is relevant to the cross section. 
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3 Ingredients for NNLL threshold resummation 



The resummation of threshold logarithms proceeds by using the factorization scale inde- 
pendence of the total cross section to derive renormalization group equations for the hard 
function Ha', which appears in eq. (11.31) . and the soft function. To define the NLL, NNLL, 
etc. approximations, we note that near threshold the usual expansion, where a s In (3 counts 
as order one, is combined with an expansion in (3, such that a s //3 also counts as one. This 
leads to a parametric representation of the expansion of the cross section in the form 

<3v = o- (0) S~] ( ^ ] exp ha. f3 g (a s ln f3) + g 1 (a s ln(3) + a s g 2 {a s In/?) + . . . 

t^V/ 9 / l " ' v ' v . ' 

(LL) (NLL) (NNLL) 

x {1 (LL,NLL); a s , (3 (NNLL); a 2 , a s f3 } (3 2 (NNNLL); . . .} , (3.1) 

which reproduces the standard structure [10] away from threshold for k = and no expan- 
sion in (3. Thus, in fixed orders, LL includes relative to the tree term all terms of the 
form 

LL a s 1 1 In 2 /?} ; a 2 s In 4 /?} ; . . . , (3.2) 

while NLL and NNLL further include all terms 

.2 fln/5 , 3 



NLL a s ln# j — ,ln d /?j> ; . . . , 

NNLL a s {l,(3x In 2 ' 1 /?}; a 2 j-^, In 2 ' 1 /3, /? x ln 4 ' 3 /?| ; . . . , (3.3) 

respectively. Note that while the LL approximation sums soft logarithms of the form 
a™ log (3 m = a s log (3 2 , . . . with n + 1 < m < 2n, this does not include all terms of this form 
at 0(a 2 ). Similarly, the NLL approximation sums soft logarithms of the form (a s log/3) n 
but 0(a 2 s log f3 2 ) terms from the interference of the one-loop hard function and leading 
soft logarithms are included only at NNLL. The NNLL terms proportional to (3 originate 
from /^-suppressed corrections to the hard functions and the soft gluon couplings. The 
former, however, vanish, since S-wave and P-wave production processes are not interfering. 
The /3-suppressed soft corrections average to zero in the total cross section at 0(a s ) (and 
probably as well in higher orders), such that no terms of the form a s (3 x In 2,1 (3 are present 
in the fixed-order expansion. 

We would now like to explain briefly the several functions appearing in the factor- 
ization formula (11.31) . We also discuss how the expansion (13 .ip is generated from this 
expression and, starting from NNLL, additional contributions of a similar factorized form. 
The derivation of the factorization formula (jl.3p in [22] relies on soft-collinear and po- 
tential non-relativistic effective field theory, and the fact that soft gluon interactions can 
be decoupled from collinear and potential fields in the leading-order effective Lagrangian 
using field redefinitions involving the Wilson lines (ll.4p and (11.51) . For the initial-state 
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partons, this redefinition is identical to that in the derivation of the factorization formula 
for the Drell-Yan process at partonic threshold [8]. The redefinitions for the final state 
particles are a generalization of those used for non-relativistic W^-bosons [23] to the case of 
a colour Coulomb force. The soft function ( 11. 6ft collects the Wilson lines arising from these 
field redefinitions. The function J# a factorizes potential effects and sums Coulomb gluon 
exchange related to the attractive or repulsive Coulomb force in the irreducible colour rep- 
resentations R a . It is defined as a correlation function of non-relativistic fields and can be 
expressed in terms of the imaginary part of the zero-distance Coulomb Green function of 
the Schrodinger equation. The hard function is defined in terms of squared short- distance 
coefficients, Hu> oc C p p,C pp ; , analogous to the corresponding treatments of heavy-particle 
pair production in e + e~ collisions [28] and of the Drell-Yan process [8]. The coefficients 
Cpp,(M, /i) encode the contribution of hard momenta to the process pp' — > HH' and are 
obtained from the scattering amplitude for the partonic subprocess evaluated directly at 
threshold [22]. 

For resummation at NLL accuracy, the required ingredients in eq. (11.31) are the hard 
function Ha/ and the soft function W^, a both at tree-level and the one-loop anomalous 
dimensions appearing in the evolution equations, with the exception of the so-called cusp 
anomalous dimension related to the leading logarithms that is required at two loops. We 
note that part of the NLL a 2 //? x In /3 term arises from the running coupling in the Coulomb 
potential and is correctly taken into account by choosing the scale to be M/3 in Jf> a . 

For resummation at NNLL accuracy the hard and soft function are needed at one- 
loop level, the cusp anomalous dimension at three loops and all remaining anomalous 
dimensions at two loops. While most of the anomalous dimensions are known to the 
required order or higher from studies of deep-inelastic scattering, the Drell-Yan process or 
Higgs production from gluon fusion, the soft anomalous dimension for pair production of 
heavy coloured particles is currently available only at one-loop [3,10,18], despite recent 
progress on massive amplitudes at the two-loop level [24,25,29]. 

In addition to the soft corrections considered in this paper, further logarithmic contribu- 
tions arise from higher-order terms in the effective Lagrangian or the production operators 
where soft gluons do not decouple after the field redefinitions. An example is the dipole 
interaction x ■ E in the potential non-relativistic QCD Lagrangian. These terms, however, 
can be treated as perturbations so that the entire expansion (13.11) can be constructed as a 
sum of terms in the factorized form x * with new hard, potential and soft 
functions, of which eq. (11.31) constitutes only the leading term in the expansion in (3. Due 
to the (3 suppression, the first correction to the leading term arises at NNLL. Therefore 
further NNLL terms in addition to (II. 3p may arise from two sources related to higher- 
dimensional terms in the expansion in (3 in the factorization formula: first, corrections 
to the Coulomb function J Ra due to subleading heavy-quark potentials contribute NNLL 
terms, which at 0(a 2 s ) are of the form a 2 log/3 [28]. Second, higher- dimensional soft func- 
tions with insertions of the x ■ E interaction potentially also contribute a 2 log 2 ' 1 (3 terms 
from the interference of a /3-suppressed soft-gluon emission with one-Coulomb exchange. 
As mentioned before, there are no corrections linear in f3 related to the hard function. 
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Finally, we mention that before the convolution with the parton distributions, the NNLL 
resummed cross section should be matched to a fixed-order two-loop calculation that is 
not yet available for the processes of interest. 

In this section we continue our investigation of the (leading) soft function and provide 
all those ingredients required for an NNLL resummation for pair production of arbitrary 
coloured particles that are related to (II. 3p except for the process-dependent one-loop hard 
functions H^, i.e. we provide the one-loop soft function W^, a and its two-loop soft anoma- 
lous dimension. In subsection 13.11 we compute the one-loop soft function for initial-state 
particles in arbitrary representations r and r' of SU(3) and a final state system in an arbi- 
trary representation R a , using the reduction of the general soft function for two final state 
Wilson lines (11.61) to the soft function (I2.23|) with a single Wilson line representing the 
final state system, achieved in section [2j In subsection 13.21 we then obtain the two-loop soft 
anomalous dimension using results of [25,29,30]. In subsection 13.31 we relate our results to 
the conventions used in the NNLL treatment of the top-quark production cross section in 
Mellin space [11]. 



3.1 One- loop soft function 

In this subsection we compute the one-loop term in the loop expansion of the soft function 

n=0 ^ ' 

The n-loop contribution to the soft function is obtained from the definition (I2.23P by 
expanding each Wilson line up to order gf 1 and keeping all contributions to the soft function 
of order g^ n . At tree level, the soft function in position and momentum space is simply 
given by 

At the one-loop order, the expansion of the Wilson lines gives rise to real and virtual 
initial-initial (ii), initial-final (if) and final-final (ff) state interference diagrams. Examples 
for diagrams contributing to the real corrections are shown in figure [2J As an example con- 
sider the initial- final state interference diagram denoted by (if) in figure [2J This arises from 
the contribution to the soft function (12.231) where the Wilson lines S^ aii (0) and S^^z) 
both contribute at order g s while all the other Wilson lines give trivial contributions. This 
diagram arises from the expectation value 

(0|T[S'^(^)]T[^ li (0)]|0>| if 

/* oo /*0 

-igsTf^ ) (iff.T^) J ds J dt (0\T[v ■ A c (v(z° + s))]T[n ■ A d (tn)}\0)\ i{ 

/0 r0 
ds / dtD + (v(z°-s)-nt). (3.6) 
■oo J — oo 
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3 3 3 

(ff) (if) (ii) 

Figure 2: Examples for diagrams contributing to the soft function 



Here we have introduced the cut gluon propagator in position space 

Dfj+ix) = (0|^(x)4(0)|0) = (-g, u )8 ab D + (x). (3.7) 

It is convenient to evaluate the integrals arising in the soft function directly in position 
space using the dimensionally regularized form of the cut-gluon propagator [31] 

= w ^7 ( 3 - 8 ) 

Atx z e [— (x + — io){x-. — ib)\ l e 

with x + = n ■ x and a;_ = n-x, and where the last expression holds for x^_ = x^ — a;_ n/*/2 — 
x + n^/2 = 0. The usual Feynman propagator in position space is given by 

Proceeding in the same way for all diagrams we obtain the one-loop soft function. 
Setting scaleless integrals to zero, only the real corrections are non- vanishing. The complete 
one-loop soft function can then be written in terms of group-theory factors C and integrals 
X as 



{aa,bf3} » 



The group theory factors are given by 



(3.10) 



Cf ] hm =2T^T { p h a 5 aP , (3.11) 

{aa,bp\ ai&i £1202 a P' v ' 

r<( if ) _ o frp(r)a r r iTi(r')a\ rp(fla)a /n io\ 



'{aa,6/3} * ^ A a 1 b 1 u a 2 b 2 T" "aibj J- a2 fc 2 y /3q 
^{aa,&j8} = ^oibi^oafcaT^™ "Tm'" = Cfl<Ai 6i<W 2 <W (3.13) 

For the initial-initial state diagrams it was used that diagrams with a soft-gluon coupling 
to two collinear particles in the same direction vanish. The integrals are given by 

l^(z Q ,a) = fi 2e (n ■ n) [°dsdtD + (z v + tn-sn) = -£fcfli e 7 B ^!£o/A 

7-oo 8tH e V 2 / 
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l {ii \z ^) = tf\v ■ n) J dsdtD + ((z - s)v - nt) = ^^r~/ B , (3.15) 

J( ff )(. ,/,) = fl 2 v 2 J°dsdt D+(v(z + s-t)) = (i^ef " ' ^ 



Here ft 2 = /z 2 e 7E / '(47r) and an infinitesimal imaginary part zq — > zo ~ ^ is kept implicit. 

It is useful to identify the amplitude A pp >{ a } with a vector in colour space denoted by 
\App>) [4,32]. More precisely the amplitude is a matrix element A pp >{ a } = ({ a }\A pp i) with 
an orthogonal basis In our notation an antiparticle in representation r is described 

as a particle in the representation f. For incoming and outgoing particles the action of a 
generator acting on particle % in our conventions is given by 



({6}|T^> = {-T^l a )A bl ... ai ... bn , incoming particle 



- ciibi 

i(Bo),« 



(3.17) 



{{b}\T°\A) = A bl ... ai ... bn Tl2 h , outgoing particle. 
Colour conservation implies the identity 

i 

In this notation the decomposition (12.131) reads 

i 

and the components of the soft function (I2.23P are expressed as W^, a = (c ( -*'- ) |W Rq \c^). 
Combining the results from eqs. (13.111) to (13.161) and expanding in e, we find that the 
(unrenormalized) one-loop soft function in the colour-operator notation reads: 

W«^(L) = - ((Tj + T 2 ) • T 3 + 2T X • T 2 ) + 2 -L + L 2 + ^) + T| (-+2L + A^j 

= (T? + T 2 ) (| + + L 2 + ^ + 2 T 2 Q + L + 2^ . (3.20) 

Here we have introduced the variable [31] 

L = 2.n(^fZl), (3.21) 

and used colour conservation ( 13.181) to arrive at the second equality. 

The components W^, a entering the factorization formula (I1.3P can be obtained from the 
above result by contracting with the elements of the colour basis cj^ Q according to (12.221) . 
As discussed in section 12.44 a large number of matrix elements are zero by construction 
and the soft matrices assume a block-diagonal form (c.f. (I2.25P ). Since the tree-level and 
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one-loop soft functions are proportional to the unit matrix in colour space and the colour 
tensors in our basis are given in terms of the Clebsch-Gordan coefficients (I2.24p . we find 
that the components ( 12.221) are diagonal due to the orthogonality of the Clebsch-Gordan 
coefficients: 

Wfr(L,fi = W«"(L,fx) 5 U , 5 RaRfj . (3.22) 

In agreement with the general results of section 12.41 these elements are non- vanishing only 
if the final state representation R a is identical to that in the pair Pj = (r a , Rp) that defines 
the basis element The diagonal elements at tree- and one-loop level are given by 

W^ )Ra (L,fi) = 1, (3.23) 

W i > V)Ra (L, /i) = (C r + C r ,) f| + -L + L 2 + ?pj + 2C Ra Q + i + (3.24) 

As shown in appendix [O the Fourier transform of this result agrees with [21] for the special 
case of the production of a colour-octet particle from gluon fusion. 



3.2 Renormalization group equations and anomalous dimensions 

In the following we will provide the evolution equations of the hard function Ha* and 
the soft function W Ra and determine the relevant anomalous dimensions at the two-loop 
level, as required for resummation at NNLL accuracy. As a result of the rewriting of the 
soft function in section [2j the soft anomalous- dimension matrix for pair production at 
threshold is identical to that of a two-to-one scattering process with two massless legs and 
one massive leg. Employing results from a recent analysis of constraints from soft-collinear 
factorization on the structure of infrared (IR) singularities of scattering amplitudes with 
massive particles [25] , we extract an analytical expression for the two- loop soft anomalous 
dimension at threshold. Explicit expressions for the resummed cross section in momentum 
space will be given in [22], but we stress that the results given here are also applicable to 
the resummation in Mellin-moment space. The precise relation to the formalism in Mellin 
space is discussed in section 13.31 

As mentioned above, in the effective field theory treatment of the factorization formula 
the hard function is defined in terms of short- distance coefficients C^, that are obtained 

from the components A®, of the scattering amplitude (12 . 13[) evaluated at threshold [22]. 
After renormalization of the ultraviolet divergences, the short- distance coefficients contain 
further IR divergences that match the ultraviolet divergences of the long-distance objects 
in the factorization formula. The IR-renormalized coefficients obtained by minimal sub- 
traction of the IR-poles satisfy an evolution equation of the form 

(M, /i) = Tij(M, ,x) Gg(M, //) (3.25) 

with an anomalous-dimension matrix whose form is constrained by soft and collinear fac- 
torization [25,33-35]. Since we have shown that the (leading) soft function is diagonal 
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for all cases relevant to hadron-collider processes, only the diagonal elements of the hard 
function, Ha = Hi enter the formula for the production cross section (11.31) . They satisfy 
the evolution equation 

-jt — Hi(M,fi) = 2ReTi(M, fi) Hi(M, jj,) (3.26) 
am/i 

with Tjj = Tj. As shown in appendix [TJ the results of [25] constrain the anomalous 
dimension to be of the form 

+ 7 f. (3.27) 

Here we have introduced the coefficient 7 CUS p by writing the cusp anomalous dimension 
for a massless parton in the representation r in the form r£ usp = C r 7 CUS p consistent with 
Casimir scaling which is appropriate at least up to three-loop order. The explicit one- 
and two-loop results for all anomalous dimensions needed in this section are collected in 
appendix IDl The coefficient 7 cusp is known to three-loop order [36], and the a s (a s rif) k 
terms are known to all orders [37]. 

Adopting the result of [25] for the structure of the anomalous dimension matrix, at 
least up to the two-loop level the anomalous dimension Yt can be written in terms of 
single-particle anomalous dimensions: 

jV =Y + Y' + Tfc- (3-28) 

The one- and two- loop anomalous- dimension coefficients Y of massless quarks, 7 9 = 7 3 , 
and gluons, 7 9 = 7 8 are given in appendix A of ref. [35]. The anomalous dimension 7^ is 
related to a massive particle in the final state representation R a in the pair P { = (r' a , R a ) 
defining the colour basis element with index i. 

It should be mentioned that eqs. (13.271) and (13.281) are derived from ref. [25], where it 
is assumed that the two heavy particles have fixed but unequal velocities, when the poles 
in e of the hard amplitude are extracted. This is different from the limit we consider here, 
where (3 — > before the limit e — > and before the loop integrations are performed, which 
corresponds to the threshold expansion as defined in [38]. The order of limits does not 
commute, and by expanding in (5 first new IR divergences appear in the hard region that 
do not correspond to UV divergences in the soft and collinear but in the potential region 
of the threshold expansion. The complete result for the scale dependence (I3.26P of the 
hard coefficient may thus contain additional terms related to the ultraviolet divergences of 
the higher-dimensional heavy-quark potentials and soft functions. In deriving eq. (I3.27P 
from the 2 — > 1 process with a single particle in representation R a we implicitly set to zero 
the scale dependence of the hard function related to the contribution from potential diver- 
gences and the higher-dimensional soft functions. This allows us to relate the anomalous 
dimensions (I3.27P and (I3.28P directly to the one of the leading soft function W^, a . 

For the case of a heavy quark an analytical result for the two-loop anomalous dimen- 
sion 7^ = 7 ff s has been extracted in [25] from the anomalous dimension of the heavy-light 
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Ti(M,fi) = -7cus P 
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2 



\C r -\- C r 



In 



AM 1 
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mC 



R c> 



quark current in SCET. In order to generalize this result to arbitrary representations R, 
we observe, following [25], that appears in the anomalous dimension of the HQET 
heavy-heavy current for a heavy particle in representation R, 

Tjn h = CWcuspCA a.) + 2 7 « s . (3.29) 

The cusp anomalous dimension 7 CU sp(/3, a) is a function of the cusp angle cosh/5 = V\ ■ t>2, 
with v\ t 2 the four-velocities of the heavy particles. The anomalous dimension V j hh for 
heavy quarks is available at two-loop order [30]. For large cusp angle the massive cusp 
anomalous dimension is related to the cusp anomalous dimension for massless particles 
according to [25] 

7cus P (/5, a s ) -> 7c US p(as) /3 + ■ ■ ■ , (3.30) 

where the remainder vanishes for /3 — > oo. It follows that the heavy particle soft anomalous 
dimension 7^ s can be obtained as one-half of the constant coefficient in the anomalous 
dimension of the heavy-heavy current f 1 3 . 2 9 1) in the limit where the cusp angle goes to 
infinity. Since the anomalous dimension of the heavy-heavy formfactor in HQET is related 
to the expectation value of a Wilson line [30] the colour structure is constrained by the 
non-abelian exponentiation theorem [39,40]. For a heavy particle in the representation R, 
at the two-loop level only the colour structures CrCa and CrTfUj appear. We therefore 
obtain the one- and two-loop anomalous dimensions for a heavy particle in an arbitrary 
representation from the result for a heavy quark by simple Casimir scaling, 

7h> = C Ra 7H,s, (3.31) 

with 7# jS = 7Q jCp. Adopting the two-loop anomalous dimension of the HQET formfac- 
tor [30] in the explicit formulation in terms of polylogarithms given in [29] we obtain 



98 2tt 2 , \ 40 „ (3-32) 



(0) 
T H ,s 


= -2, 


-Y (1) 


= -c A 



9 3 + 4 ^ + ¥ r pn/ , 

where the loop expansion of the anomalous dimensions is defined as in (ID. 101) . The one- loop 
expression in (13.321) agrees with the well-known result that the one-loop soft anomalous 
dimension is proportional to the quadratic Casimir of the final state system [3,10,18]. The 
two-loop expression for the anomalous dimension (13.281) for heavy-particle pair production 
in an arbitrary colour representation is a new result. The Casimir scaling of the single- 
particle anomalous dimension (13.311) was also noted in the published version of [25]. 

The factorization scale independence of the hadronic cross section in the threshold 
region can be used to obtain the evolution equation of the soft function as [22] 



d 



-W^(L) = ((r^ usp + r^ sp )L - 2 7 fe) W*~(L), (3.33) 



dln/j, 

where the anomalous dimension of the soft function is obtained from the anomalous di- 
mension of the hard function by adding anomalous dimensions 7^ ,r entering the evolution 
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equations of the parton distribution function for a parton in the representation r in the 
x — > 1 limit (see ID.161) : 

7^ = 7^ + 7^ + 7^'- (3-34) 

The anomalous dimensions r y^' r for quarks and gluons are available up to three-loop or- 
der [36]. The one- and two- loop results are collected in appendix [Ql We note that eq. (I3.34p 
is derived under the assumption that the Coulomb function J Ra is scale- independent, 
dJR a /d\n fi = 0, which is no longer true starting from NNLL in the non-relativistic expan- 
sion. However, as remarked above this type of scale dependence is related to terms that 
we consistently dropped in the derivation of eqs. (13.271) and (13.281) . so the two-loop soft 
anomalous dimension of the Wilson line considered here is unaffected by these complica- 
tions. 

Analogously to the anomalous dimension j[ of the hard function (13.281) . at least up to 
the two-loop level the anomalous dimension of the soft function (I3.34[) can be written in 
terms of separate single-particle contributions 

lw7 = 72> + 7l + < (3.35) 

with 

Y s =Y + 7*+ (3-36) 

Up to the two-loop level, the anomalous dimensions for quarks 7f = 7^ and gluons 7f = 7® 
are related by Casimir scaling, 

7 I = C rls , (3.37) 
where the one- and two-loop coefficients are [8]: 

7< 0) = 0, 

m „ ( 404 IItt 2 \ m fll2 2tt 2 \ (3.38) 

7s — Ca y ^ ig~ + l^CaJ + T F Uf y— — j . 

The evolution equations (13.261) and (13.331) generalize the corresponding equations for the 
Drell-Yan process [8] and Higgs production [41] to processes with a heavy particle pair in 
the final state. The evolution equation for the soft function in momentum space involves 
distributions and can be solved in Mellin-moment space [31] or directly in momentum 
space [6-8] using a Laplace transform. Using the two-loop result ( 13.321) . the anomalous 
dimension in the evolution equation (13.33!) i s known with the accuracy required for NNLL 
resummation. 



3.3 Relation to the formalism in Mellin-moment space 

In the applications of resummation in Mellin space [9-13, 18, 19] the Mellin moments of the 
partonic cross section with respect to the variable p = 4M 2 /s are written in the form 

■N fn/r2 ,.\ — / j„JV-1* (a 7iyf2 



= ^ N R SM\p)g^ Ra {M\p) exp (G^(M 2 , /, 2 )) . 
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Here & pp / j R a is the partonic cross section for the production of a heavy-particle final state 
pair in the representation R a , u^, N Ra are the Mellin moments of the Born cross section, the 
matching functions g p ^ collect the iV-independent corrections, and the exponent G pplRa 
has the form 



G^, Ra (M 2 ,^) = / dz- / % ((Ap (a s (q 2 )) + A p , {a s {q 2 ))) 

Jo 1 — z Jp? Q 



AM*{l-z) 2 j 2 



+ D«'-*hh> K(4M 2 (1-^) 2 )) 



(3.40) 



The coefficients A p contain the effect of collinear radiation off the incoming partons and are 
identical to the cusp anomalous dimension T r cusp . The coefficient D describes soft radiation 
and can be written as a sum of terms related to the incoming partons and the final state 
system 

D^hh'M = \ (D p (a s ) + D p ,(a s )) + D^ H ,(a s ). (3.41) 

Applications of resummation in Mellin space often use an equivalent form of (13.391) where 
iV-independent terms contained in the large-iV expansion of (13.401) are not included in 
the exponent. This implies a redefinition of the matching functions gp pRa (M 2 ,fi) but no 
change in the A and D coefficients. 

The coefficient D pp a ,_^ HH , is not identical to the anomalous dimension of the soft function 
7^* because the Mellin transform of the resummed cross section in momentum space 
involves the fixed-order soft function at the scale fi s ~ M/N [8,27] while in the Mellin- 
space formula (13.391) all the iV-dependent soft corrections are exponentiated. Furthermore 
the form of the exponent is different in the two approaches. The relation of different 
forms of exponentiated expressions in Mellin space is discussed in [7,8,37,42,43]. But for 
heavy-particle pair production near threshold, the structure of the resummed expressions 
in Mellin space and momentum space is identical to that for the Drell-Yan process, where 
the relation between the two formalisms was obtained in eq. (71) of [8]. This then implies 
a relation of the coefficient (I3.4ip to the anomalous dimension of the soft function 
given by 

e 2 7B v r(1 + 2V)D^ HH , =2 7 ^ + 2Vlnsf«(0,/i) 

e 27 E v r (! + 2V) - 1 
V 

where 

d (3(a s ) d 



dliafi 2 2 da s ^ 

We also introduced the Laplace-transform of the soft function [6, 7] with respect to the 
variable s = l/(e 7B /i e P//2 ) 

POO 

g?»(p ) Ai)=/ due-"* Wf>,/i), (3.44) 
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where we have defined the MS-renormalized soft function W i a . Since the soft func- 
tion in position space W(zo,f-i) depends on the arguments solely through the variable 
/2) = e L l 2 with L defined in fl3T2T|) . it is easy to see that the function s(p) is 
obtained by simply replacing L — > —p in the MS-renormalized result for the soft function 
in position space. Expanding the relation (13.421) counting V ~ a s we obtain the terms 
relevant to determine D at the two-loop level: 

D^hw = 27gT, + 2Vlnaf-(0, ai) - y V(I^ usp + T^ sp ) + 0(a 3 s ). (3.45) 

Using the loop expansion of the soft function (13.41) . the anomalous dimensions (ID. 1Q|) 
and the beta function, 



n+l 

dfi 4vr 
with p Q = 11/3 Ca — 4:/3Tpnf, we find up to the two-loop level 



(3(a s ) = ^ = -2a s £ ) & (3-46) 



D { X a HH ' = 2 ^wT a > ( 3 - 47 ) 

d^ hhi = 27^ - 2/?o (f"-{o,ri - y(rS + rg)£)) 

= 2 7 «f Q + ft (vr 2 (C r + C r ') - 8C R J , (3.48) 

where we have used Tculp = 4C r and the one-loop soft function ( 13.241) in the last step (with 
L = and the 1/e poles discarded due to the MS subtraction). 

Using the decomposition (13.35)) of the anomalous dimension of the soft function and 
applying (13.481) to Drell-Yan and Higgs production, in which case r = r', Cji a = (colour 
singlet final state) and r ^\ Vi = 2 7 J due to (I3.35p . we obtain the coefficients D q and D g 
introduced in (I3.4ip : 

D p = (4 7 ^ + 2vr 2 a/3o) (^) 2 + C{al). (3.49) 



The result for D q agrees with eq. (72) of [8] and the explicit expression obtained using (£5. 
agrees with the one used in eq. (A. 3) of [11]. 

Turning to the soft anomalous dimension related to the heavy particle pair, we find 
that the relation of the coefficients appearing in the Mellin-space approach to the soft 
anomalous dimension obtained in (13.321) is given by 

Df*« = 27^, (3.50) 
0^ = 241^-8^0^. (3.51) 

The one-loop coefficient D^^T = —ACn a is in agreement with previous results for top- 
quark and gluino production [3, 10, 18]. The result for the two-loop coefficient is new 
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and shows that at NNLL level there is a non-trivial relation between the soft anomalous 
dimension and the resummation coefficients if eq. (13.391) is used for the resummed cross 
section. Our result for the two-loop soft anomalous dimension (13.321) allows to obtain 
for the first time the two-loop coefficient for soft radiation off a massive particle pair at 
threshold for all possible colour states of the heavy-particle system: 

D gg. = -cs, (x - £ + 8C3 ) + f c ^ n >- < 3 - 52 > 

Our result (I3.52p differs from the result quoted in [11] for the special case of a colour-octet 
final state by a term 8Ca [Ca(1 — Cs) ~~ A)]- However, the result in [11] was obtained from 
an analysis of the singularities of the two- loop massive quark form factor [44,45] in the limit 
of light quark masses using the unjustified - and as it turns out incorrect - assumption 
that the factor of proportionality of the two-loop and one-loop soft anomalous dimension 
matrices in the massless limit is identical to the one at threshold. 

4 Conclusions and outlook 

We have performed a detailed study of the soft function relevant to threshold resummation 
of production processes of heavy coloured particle pairs at hadron colliders. We have given 
a precise formulation of the physical picture of soft-gluon radiation coupling to the total 
colour charge of the heavy-particle pair. This has allowed us to construct a colour basis 
that diagonalizes the soft function to all orders in perturbation theory. Explicit expressions 
for all production processes of top quarks, squarks and gluinos have been provided. We 
have calculated the one-loop soft function for arbitrary colour representations of initial and 
final state particles and used recent new insights into soft-collinear factorization to obtain 
the two-loop soft anomalous dimension. This supplies the process-independent ingredients 
for NNLL resummation of threshold logarithms in arbitrary production processes of heavy 
coloured particles at hadron colliders. A complete NNLL resummation in the sense of 
eq. (13.1 p further needs the colour-separated one-loop short-distance coefficients, as well as 
the summation of logarithms associated with subleading terms in f3 in the non-relativistic 
expansion. 

In a subsequent publication [22] we will give a derivation of the factorization for- 
mula (ll.3p that demonstrates the factorization of Coulomb gluon exchange from soft-gluon 
radiation using field redefinitions in an effective field theory. The formula thus provides a 
theoretically clearly defined separation of hard, Coulomb and soft effects and can be used 
for a combined resummation of soft and Coulomb gluons. This will be discussed for the 
case of squark-antisquark production, where the effect of Coulomb resummation may be 
of similar order as the effect of soft gluon radiation [19]. 
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During the review process of this paper refs. [50] appeared, in which the authors reported 
a non-vanishing result for the three-particle contributions to the two-loop soft anomalous 
dimension of amplitudes with two massive particles and any number of massless particles 
near threshold, which is also not diagonal in the colour basis discussed in the present paper. 
We emphasize that these findings are not in contradiction with the results reported here. 
The authors of refs. [50] calculate the 1/e poles away from threshold and then take the 
limit (3 — > 0, resulting in logarithms of (3 in the anomalous dimension and non-zero three- 
particle correlations. In the present approach, where the expansion in (3 is constructed 
directly within the non-relativistic effective field theory framework, the extra logarithms 
are related to the potential region and the three-particle correlations to higher- dimensional 
soft functions, both belonging to the possible NNLL terms mentioned in the text that arise 
from subleading heavy-quark potentials, and from 0{f3) terms interfering with Coulomb 
singularities. Thus, the derivation of the diagonal colour basis and two-loop anomalous 
dimension for the soft function with two equal heavy particle velocities discussed in the 
present paper remains valid in the light of the results of refs. [50], but should not be 
expected to apply to higher-dimensional soft functions. (In fact, the emission of a soft 
gluon due to the subleading x ■ E interaction mentioned in section [3] leads to a change 
in the colour state of the heavy-particle pair, and the corresponding soft function must 
therefore be off- diagonal.) Finally, let us mention that while the results of refs. [50] may 
be used to determine all logarithmic terms at 0{a 2 s ), the summation of NNLL logarithms 
in higher orders in the strong coupling expansion requires an analysis of soft and potential 
divergences directly at threshold. 

A Technical details on the colour structure 
A.l Construction of the colour basis 

To see that the basis can always be chosen as in (12.171) . we use the completeness rela- 
tion (12.31) of the Clebsch-Gordan coefficients to decompose the scattering amplitude (I2.13P 
according to 



Note added 




(A.l) 



r a ,Rg 
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with coefficients 

4# = cXbAAm- ( A -2) 

The fact that the amplitude A{ a y is colour conserving, i.e. satisfies an identity analogous 
to (12.151) . and the invariance condition of the Clebsch-Gordan coefficients imply that the 
coefficients ( 1A.2j) are invariant tensors under transformations in r a (g) Rp: 



A aP = U^A, 5 Ufc\ (A.3) 

Hence according to Schur's Lemma A a p is non-vanishing only if r a and Rp are equivalent 
irreducible representations, i.e. belong to one of the pairs Pj = (r a , Rp) of equivalent 
representations appearing in the decompositions of the initial and final state system into 
irreducible representations introduced above eq. (12.161) . In this case A a p = A^ 1 ' S a p for 
some coefficient Ay> . We have therefore derived the decomposition 

Pi 

of the amplitude. Using the orthogonality of the Clebsch-Gordan coefficients we find that 
the constant A^ is, up to normalization, given by the components of the amplitude in the 
basis f l2~TTj) 

The identity ( 1A.4I) is therefore precisely the decomposition of the amplitude into the ba- 
sis (12.171) for the coefficients. 

A. 2 Identities for Wilson lines 

In order to rewrite the soft function in terms of the Wilson line for a single heavy particle 
in a representation R a we used the identity (12.191) and its complex conjugate 

c-(R')t o(-R)t r R a* — r R «* c( fia )t t a (\\ 

J v,b 2 a 2 J v,b 1 a 1 Ly aa 1 a 2 ~ ° f3b 1 b 2 °v,f3a ' l^- u J 

To prove identity (12.191) we observe that the Wilson line (jl.5p solves the differential 
equation 

(v-D)Si R) (x ) = (v-d-ig s vA a T a ^)Sl R \x ) = 0. (A.7) 
Then using the relation 

riRa frpa(R)r r rp^R'A _ rpa(R a ) ^-y_R a /. o\ 

^aa ia2 ^""-aibi °a2f>2 ' u a 1 b 1 1 a 2 b 2 J ~ 1 a/3 °/36if>2' V 1 ^- ) 

obtained from the invariance condition (12. 7p for infinitesimal transformations, we see that 
the left-hand side of (I2.19P also satisfies ( 1A.7h for a Wilson line in the representation R a , 
if the Wilson lines Si, and Si, satisfy the analogous definitions in their representations: 

(11 n\ r* Ra q( r ^ c( r/ ) 

yu ■ U )a/3^p ai a 2 v,a 1 b 1 °v,a 2 b 2 
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(~<Ra 



((v ■ D)Sl R \ lbl Sif 2b2 +Si R l bl ((v ■ D)Si R \ 2b2 = 0. (A.9) 



Since the Wilson line ( 11.51) is the unique solution of the differential equation flA.71) with 
boundary condition lim^oo S v (sv ) = 1, eq. ( 12. 19f) follows. 

B Clebsch-Gordan coefficients, colour bases and pro- 
jectors for gluino and squark production 

B.l Squark-antisquark (top-antitop) production 

We collect here the Clebsch-Gordan coefficients, projection operators and basis tensors for 
quark- ant iquark and gluon-gluon initiated production of a 3(g) 3 final state. The projectors 
on the singlet and octet final state representations are independent of the production 
channel and have been given already in (12.121) : 

P (1) = — 5 6 

r {a} N ^a ia2 O a:iail 

7}(8) O rpa rpOL 

{a} ~ a 1 a 2 1 a 4 a 3 - 

The indices take the values a; L £ {1, 2, 3}. 
B.l.l Quark- antiquark fusion channel 

For completeness we repeat the results for the quark-antiquark channel given already 
in (12. lip and (12.180 . The Clebsch-Gordan coefficients for the two representations in the 
decomposition 3 ® 3 = 1 + 8 are 

aia2 (B.2) 
r-<(8) _ x /oT a 

^aaia,2 v 0201 

with a G {1, . . . ,8}. The basis elements for the colour structure of the hard production 
process corresponding to the two possible combinations Pi = {(1, 1), (8, 8)} are given by 

(1) _ J_r c 

{a} ~ at aia2°a 3 a4 > 

o (B.3) 

(2) = Z rp/3 rpP 
{a} ■\fD~A "314 

with D A = iV 2 - 1 and N c = 3. 
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B.1.2 Gluon fusion channel 

For the production of a 3 <8> 3 final state from gluon fusion there are three possible combi- 
nations of equivalent initial and final state representations: 

Pi e {(1,1), (8 5 ,8), (8 A> 8)}. (B.4) 

The Clebsch-Gordan coefficients for combining two particles in the adjoint into a singlet, 
a symmetric and an antisymmetric octet are 

^(i) _ 1 x 
CK 8s ) = - D a CR 51 

(](&a) = 1 ™ 

where all indices run from 1 to 8. Here we have defined = if aiaa2 in terms of 

the SU(3) structure constants, the symmetric invariant tensor -D" 2 = d aaia2 , and the 
coefficient £?p = = A appearing in the relation Trf/}"/}^] = 4B F 5 a /3. From the 

definition (12.171) we obtain the same colour basis that has been found to diagonalize the 
one-loop soft anomalous dimension matrix [3] 



'A 



{a} y/N c D/ 2 4 ' 



{a} y/2D A B F a2ai a3a4 ' ( bj 

r ( 3 ) — , / _ F Q T« 

{a} y jy (12!ll «3«4' 

B.2 Squark-squark production 

For quark-quark initiated processes qq —>■ qq the initial- and final-state systems are either 
in the 3 or 6 representation of SU(3) since 3 ® 3 = 3 + 6. We denote the 6 by a symmetric 
double index a = (ccia^) with oti G {1, 2, 3}. The Clebsch-Gordan coefficients are given by 

f7(3) _ 1 , 

V2 (B.7) 



where the indices and the index a for the case of the 3 can take the values 1 to 3. For 
the sextet representation the normalization condition uses a symmetrized definition of the 
Kronecker-delta for the double indices: 

c(6) 1 



Kf3 — n(<Wl<W2 + <Wa<Wl)- (B-8) 
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The projectors on the two representations needed for the decomposition of the Coulomb 
Green function are given by 

Since the tensor product 3 ® 3 is of two identical representations and since the Clebsch- 
Gordan coefficients are real, the elements of the colour basis coincide with the projectors, 
up to normalization: 

VWc-l) (B1Q) 

c { a } = ^J 2N (N + 1) ^ aia3 ^ Q2a4 ~^ ^ a ^ a ^ a ^ a a) ■ 
B.3 Gluino-squark production 

For gluino-squark production qg — > qg the relevant representations appear in the decom- 
position 3®8 = 3 + 6 + 15. The Clebsch-Gordan coefficients are given by (for N c = 3) 

Z"f(3) _ _ rpQ2 
^ ' aa\a 2 — L aa\ i 

✓7(15) _ _ ( x rpa 2 1 r ma 2 _ ma 2 _ ^ r ma 2 ] 

w «aia 2 ^2 I " a l a i a 2 a 3 ~ "a 2 ai ->- aia 3 ^ "0103 - 1 - a 2 a\ ^ u a 2 az ± a 1 a 1 J 1 

where ai £ {1, 2, 3} and <22 £ {1, . . . , 8}. For the 6 we use the same double-index convention 
as in (IB. 71) . For the 15 we have introduced a triple index a = (aia^a^) where the first 
two indices transform in the 3 and the last index transforms in the 3 representation. The 
Clebsch-Gordan coefficient is symmetric under the exchange ot\ <-> «2 and vanishes upon 
contracting a.\^ with a 3 . In the normalization of the coefficient C^ 15 -* we use a Kronecker 
delta that has the same symmetries as the coefficient in both index triples: 

111 

&a0 = 2 ScufotiaaPa ~ g ^aias^/fe'Ws ~ g^aias ^aa/Ji <Ws + ( a l ^ (B.12) 

The projectors can be written as (for iV c = 3) 

p( 3 ) _J (rpa 2 rpa4\ 

iffy = l^ ia3 S a2a4 - l ~{T a *T a % iaz - (T a *T a % ia3 , (B.13) 

p( 15 ) _£ £ _ _ (rpa 2 rpa4\ _|_ (rpa,4rpa 2 \ 

r {a} ~ 2 aia3 a2CL4 /aias r l J J ,Mia3- 
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The colour basis that diagonalizes the one-loop soft function is related to the projectors 
by complex conjugation and a normalization factor 1/ A/dim(r a ): 

' ( r pcLA r pa2\ 

\ A A /cyan 



.(1) 
'{•} 

.(2) 
>} 

.(3) 



(^ a i a 3 ^a 2 a 4 



1 



1 

1 r ) V2 



(T a4 T a2 ) a3ai - 2(T a2 T' 



a 4 \ 



10301 J j 



«4 n 



(B.14) 



B.4 Gluino pair production 

The projectors on the several final-state representations for gluino pairs appearing in the 
decomposition of 8 (g> 8 are the same for the quark-antiquark and gluon-induced processes 
and can be obtained from [46,47]: 



p(l) 

W 


- h A 

g°aia2°a3a4; 






p(8s) 
>} 


_ _ D a jja 

g -^0402-^0403 5 






P (Sa) 
W 


— _ P Q F a 

^ 0102 04 03 ' 






p(10) 
W 


- 1 ^ A 

^ 1 0,10,3° 1120,4 


^0104^0204 


g 111 (12 0403 


p(10) 

>} 


1 /a A 


^0104^0204 


^ pa pa 

g 0\02 0403 


p(27) 
>} 


- 1 fx A 

2 I °ai 03° 0,204 


"I - ^aia 4 "a2a 3 


X A A 

^°aia 2 °a 3 a4 



+ D a pa + pa pa 

1 0301 0402 1 0301 0403 



0301 0402 0301 0402 i ' 



—D a D a 

^^020^0304 



(B.15) 



B.4.1 Quark-antiquark fusion channel 

For quark-antiquark induced processes, the pairs of equivalent combinations of initial and 
final state representations are 

Pi G {(1,1), (8,85), (8,8 A )}. (B.16) 

The Clebsch-Gordan coefficients that combine the initial state quarks into a singlet and an 
octet have been given in flB.2j) while the coefficients for the final state appeared in (IB. 51) . 
The resulting colour basis is related to ( IB. 61) by exchanging initial and final states: 



c 



>} 

.(2) 
'{0} 

(3) 
{«} 



1 



0304 1 



^ rpQ> T) a 

^/2D aBf a2CLl 3(14 ' 



(B.17) 



N C D A a2ai a3a4 ' 
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B.4.2 Gluon fusion channel 

For the production of two gluinos from gluon fusion, the allowed pairs of initial and final 
state representations are given by (12.161) : 

Pi G {(1,1), (85,85), (8 A ,8 5 ), (8^,8^), (8 S , 8 A ), (10,10), (10,10), (27,27)}. (B.18) 

The Clebsch-Gordan coefficients for the singlet and octet representations are the same as 
in (1B.5H . Since according to (12.171) the colour basis elements for the production of a 10, 10 
and 27 are, up to normalization, the complex conjugate of the projectors given in (IB. 151) . 
we do not need the lengthy Clebsch-Gordan coefficients for these representations. The 
basis consists of an operator corresponding to the production of a singlet: 

C \a\ = 7T - ^ai 02^03(14; (B.19) 

four operators corresponding to the different combinations of 85 and 8^, 



{o} in /TT~ 0201-^0304) 



1 



Q 



(3) __ _ pa p,a 

{a} 2^/B F N c D A a2ai aaa4 ' 



(4) __ ' pa pa 

{a} at /TT~ 0201 1 aza^i 



1 



(B.20) 



(5) _ pa pa 

{a} 2^/B F N c D A a2ai a3a4 ' 



two operators corresponding to the 10 and 10 

.(6/7) _ 1 



4V10 

and one operator for the 27: 



°a\a-i°aia^ a\a^fia2a^ a 2 ai^ a 3 a 4 - 1 - (^0301^0402 ' 0301 0402/ 



(B.21) 



C {a} Qy^ ^ a i a 3^ a 204 + ^0104^0203 ^^0102^0304 a 2 a 1 ^ a^a^j • (B.22) 

C Fourier transform of the soft function 

The Fourier transform of the soft function (11.7)) enters the factorization formula (11.31) while 
in section 13.11 we calculated the one-loop soft function in position space. The momentum 
space result can be obtained by inserting the Fourier transforms of the basis integrals 

r(- e ) i _ i /cj x 



^"■"^-s^^U 1 (ai) 



30 



X (ff V,/i) 



n 



8vr 2 rf-2e) 



1 (u 



(l-2e)o; V/x 



-2c 



and X^ fS) (to, ji) = —1/2 T^ n \oj, fi) into the result (|3.10p . For unstable heavy particles, the u 
integral in (jl.3p extends to infinity and the integrals have to be expanded in e in the sense 
of modified plus-distributions [48]. The Fourier transform has been calculated using 



°° teo e iu>zo/2 ( ^f0A* 
-oo 4vr 



r(— q) w \/i 



B(u), 



(C.2) 



where we recall the prescription zq ^ Zq — i5. 

To compare the one-loop soft function (I3.24p to the result obtained in [21] for the 
special case of the production of a colour octet scalar from gluon fusion, we compute the 
Fourier transform as a function of uj = M(l — z) and expand the basis integrals into 
plus-distributions, as appropriate if the decay width is neglected as in [21]: 



X^(M(l-z)^) 



8n 2 M 



21n 2 fA)-^ 
MJ 4 



1 




"ln(l - z) 




1 - z_ 


1-z 


+- 



l^(M(l-z)^) 



1 



8vr 2 M 

Here we have used the identity 



5(1 -z) 



i + 21n(A 



+ 2 -2 



'1-z) 



(C.3) 



9(1 -z) 



,-l-2e 



1-2 



- 2e 



ln(l - z) 

1 - z 



(C.4) 



for distributions on the interval [0, 1]. 

The soft function in momentum space obtained from (13.241) is therefore given by 



(l)-Ro 



(M(l-z)) 



2 

M 



+ C 



(C r -f- C r 

1 



6(1 



1 - z 



+ 4 



ln(l - z) 



M 



71" 
T 



Ra 



6(1 -z) 



1 - z 
2 1- 2 



1 - z 



(C.5) 



This reproduces eq. (40) in [21] by setting C r = C r i = Cn a = Ca = N c and multiplying 
with a prefactor a s M/(A7r) to account for our definition of W^ Ro> as coefficient of a s / (Ait) 
and the different normalization of the leading-order soft function: Wf )Ra (M(l - z)) = 



6(M(1 - z)) = 1/M6(1 - z), whereas sf /P (M(l - z)) = 6(1 - z) in [21]. 



s/p\ 
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D Anomalous dimensions 



In ID. II we derive the anomalous dimension ( 13. 27ft of the hard function. The explicit one- 
and two-loop results for the anomalous dimensions are collected in ID. 21 

D.l Anomalous-dimension matrix of the hard function 

From [25] we find that the UV regularized, minimally subtracted short-distance coefficients 
for a general 2 — ► n scattering process including massless and massive partons obey an 
evolution equation 

\C({k}, {m}, /,)) = T({A;}, {m}, /i) \C({k}, {m}, /i)) , (D.l) 



cHn/i 



where {k} = {ki, . . . , k n } and {m} = {mi, . . . ,m n } denote momenta and masses of the 
scattered particles and we use the colour-state formalism [32]. For general scattering 
processes the colour structure of the matrix T can be quite complicated with two- and 
three-parton colour correlations contributing at two-loop level. 
The two-parton correlations take the form (c.f. eq.(10) of [25]) 



r(W,{mU)|^«^^rT»P ^ ) +E^+E^fc 

~ E WW + E T ' ■ T ^ ^usp In (^) . (D.2) 

(J, J) 1,3 ^ 

Here indices i (I) denote massless (massive) partons and the notation indicates un- 
ordered tuples of distinct parton indices. We also defined = 2aijPi-pj + iO with = +1 
if partons i and j are both incoming or outgoing and = — 1 otherwise. 

The anomalous dimensions related to light partons are collected in appendix ID. 21 The 
heavy-particle soft anomalous dimension jjj s for arbitrary SU(3) representations is given 
in (13.321) . The cusp anomalous dimension 7 CU sp(/3j-j) in the terms involving two massive 
partons is a function of the cusp angle j3u = arccosh(— su/2mimj) and can be obtained 
from the heavy- heavy quark formfactor in HQET (13. 29ft . For large cusp angles the cusp 



anomalous dimension satisfies (I3.30f) while for small cusp angle as in pair production at 
threshold [25] 

7cus P (/3) ^ -2 7 Q /C F = 1h, s - (D.3) 

For processes with two heavy final state particles three-parton correlations involve the 
colour structures 

f abc T°T b jT c K , (D.4) 

while a two-loop analysis [24] and soft-collinear factorization [25] show that three-parton 
correlations with two light partons and one heavy parton proportional to f abc T^T^T^ are 
absent. 
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In the following we evaluate (ID. 21) for a 2 — > 1 process with a massive final state 
particle in an arbitrary SU(3) representation, corresponding to the reduction of the lead- 
ing soft function obtained in section [3 As discussed in section [3] this corresponds to 
dropping contributions to the anomalous dimension related to higher-order potentials and 
higher-dimensional soft functions. For the three-parton process with a single heavy particle 
the structure (ID. 41) cannot appear. Furthermore, for a two-to-one process the structure 
/ a T"T^T^ vanishes by colour conservation alone. We therefore conclude that in or- 
der to obtain the two-loop anomalous dimension of the hard function in the factorization 
formula fll.3ft corresponding to the (leading) soft function W u , a in the same equation, it 
is sufficient to consider the two-parton correlations ( ID. 21) . since all higher multi-particle 
correlations do not contribute. 

The result (13.271) for the anomalous dimension of the hard coefficient in heavy particle 
pair production at threshold can be obtained from the general expression (ID. 21) using the 
reduction to a two-to-one process with a single particle in the representation R a with mass 
2M = (ma + ruff/) and momentum P = p\ + p 2 , derived for the soft function in sec- 
tion [2j (This reduction can be performed analogously for the soft function (0\S n S n SlSl\0) 
appearing in the factorization of the amplitude to which eq. ( ID. 21) applies.) At threshold 
ki ■ k 2 = P ■ ki = 2M 2 and the two-loop anomalous-dimension matrix becomes 



+ £T,.T,Wn(^y 



(D.5) 



7cusp 



Tl ■ T 2 In ( i- 1 - To In (-^- 

V 4M 2 + i0y 3 V2M 



Ra 



+ 7' +7' +7^: 



where we have used colour conservation Ti + T 2 = — T3. Further using the identity 

2Ti-T 2 = T^-Tj-T^ (D.6) 
leads to the result for the anomalous dimension quoted in (13.271) : 



T({k,P},M,ii) = ^ lcusp 



{C r -\- C r 



In 



4M 2 
1^ 



/7T 



mC 



Ra 



+7 r +7 r '+7£- (D.7) 



It is interesting to see how the same result is obtained by applying the general ex- 
pression (ID. 21) directly to the four-particle process at threshold, omitting three parton 
correlations as discussed above. Using the fact that the heavy-particle cusp anomalous 
dimension near threshold simplifies according to (ID.3j) and pj ■ ki = mj\^/2 m mjM one 
finds 

,2 



T({k,p}, {m}, / u)|2_p ar ton = Ti • T 2 7 cusp In (~^r 



ko + iO 



+ 7 r + 7 r ' +7L + 7£ 



+ 2T 3 -T 47 ^+ Ti-T J+2Tcusp ln(^T:) 



33 



= Ti • T 2 7cus P In + l r + / + 2 T 3 • T 4 lHjS 

+ (T x + T 2 ) • (T 3 + T 4 ) 7 cusp In (^) + 7 £, s + 7?.- ( D -8) 

Using colour conservation and the analog of ( ID. 61) for a four-particle process, this expression 
simplifies to 

T({k,p}, {m},/i)|2_ par ton = ^7cusp 

+ 7 r + l r ' + [(T 3 + T 4 ) 2 — T 2 — T4] 7 h, s + ll s + 7^-(D.9) 

Since for a final state pair in a representation (T3 + T4) 2 = Cn a , we obtain the 
same result as in the three-point calculation, provided the heavy-particle soft anomalous 
dimension satisfies Casimir scaling, 7^° = C^^f . This gives a second argument for 
Casimir scaling of the soft anomalous dimension, in addition to the derivation from the 
HQET formfactor given in the main text. 



(T 2 + T 



In 



4ikP 
1? 



ITX 



+ z7rfT 3 + T 4 



D.2 Explicit results for the anomalous dimensions 

In this appendix we collect explicit results for the one- and two-loop anomalous dimen- 
sions that are already available in the literature. We define the expansion of the various 
anomalous dimensions in the strong coupling constant by 



7 



a n^ 1 



;d.io) 



The explicit one- and two-loop results for the cusp anomalous dimension are given by 

(D.ll) 



7 (°) = 4 

'cusp ' 



7 W =4 

/ cusp 



Y67 









20, 
9 



with T F 5 ab = Tr [T a T b ] and T F = 1/2. 

The one- and two-loop anomalous-dimension coefficients 7** of massless quarks, 7 9 = 7 3 
and gluons, 7 9 = 7 s , are given by [35] 

<<°)« = - 3CV, 



7 (1)9 = Cl (-- + 2vr 2 - 24C 3 ) + C A C F {- 

/130 2tt 2 \ 
+ C F T F n f (- + -y 

7 (% =-f3 = -YCA + ~T F n f , 



961 llvr 2 



54 



6 



+ 26C; 



(D.12) 

(D.13) 
(D.14) 



34 



= (-f + + 2 C3 ) + c.r^ (f - £) + 4C F r,», (D, 5 ) 

The anomalous dimensions 7^ ,r appearing in the anomalous dimension of the soft func- 
tion (I3.34p are defined by the evolution equation of the parton distribution function for a 
parton p in the representation r in the x — > 1 limit, 

f/ ■/ p/ ,(^) = 2/ r (a s )/ # (^)+2r (a s ) - ^ + ... . (D.16) 



Results for quarks, 7*' 3 = 7*, and gluons, 7*' 8 = 7 s , are available up to the three-loop 
order [36]. The explicit values in the notation used here are given by [8,41]: 

7 (W = 3CV, (D.17) 

7 W = 2vr 2 + 24 C3 ) + C A C F + ™f - 12(a) 

-c^/fi+^y (D.18) 



3 9 
3 ~ 3 



^ 0)B = Po = ^C A -^T F n f , (D.19) 



1 {1)B = 4Cl (| + 3Cs) - ™C A T F n f - 4CW>n,. (D.20) 
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